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Abstract 

We consider different aspects of the problem of cosmological singularity 
such as the BKL oscillatory approach to the singularity, the new features 
of the cosmological dynamics in the neighbourhood of the singularity in 
multidimensional and superstring cosmological models and their connec- 
tions with such a modern branch of mathematics as infinite-dimensional 
Lie algebras. The chaoticity of the oscillatory approach to the cosmolog- 
ical singularity is also discussed. The Conclusion contains some thoughts 
about the past and the future of the Universe in the light of the oscillatory 
approach to the Big Bang and the Big Crunch cosmological singularities. 

1 Introduction 

Many years ago, in conversations with his students, Lev Davidovich Landau 
used to say that there were three problems most important for the theoretical 
physics: the problem of cosmological singularity, the problem of phase transi- 
tions and the problem of superconductivity pQ. Now we know that the great 
breakthrough was achieved in the explanation of phenomena of superconduc- 
tivity [2] and phase transitions :3]. The problem of cosmological singularity 
was widely studied during the last 50 years and many important results were 
obtained, but it still conserves some intriguing aspects. Moreover some quite 
unexpecting facets of the problem of cosmological singularity were discovered. 
Isaak Marcovich Khalatnikov, who was one of the pupils of L.D. Landau, had 
brought a significant contribution into the discovery and elaboration of different 
aspects of the problem of the cosmological singularity and the chaos, arising in 
the process of the asymptotic approach to this singularity. 

In our review, published 10 years ago [4 in the issue of this journal, dedicated 
to the 90th anniversary of L D Landau birth, we had discussed some questions 
connected with the problem of singularity in cosmology. In the paper dedicated 
to 100th anniversary of L D Landau [5] we have dwelled on relations between 
the well-known old results of these studies and the new developments in this 
area. 

In the present paper, dedicated to 90th anniversary of I.M. Khalatnikov, I 
shall give a brief review of some old and new ideas, connected with the develop- 



1 



ment of the theory of the asymptotic approach to the cosmological singularity, 
and shall try to argue why it could be of interest not only for physicists and 
mathematicians, but, perhaps, for even more wide audience. 

To begin with, let us remember that R Penrose and S Hawking [5] proved 
the impossibility of indefinite continuation of geodesies under certain conditions. 
This was interpreted as pointing to the existence of a sigularity in the general 
solution of the Einstein equations. These theorems, however, did not allow for 
finding the particular analytical structure of the singularity. The analytical be- 
haviour of the general solutions of the Einstein equations in the neighborhood 
of singularity was investigated in the papers by E M Lifshitz and I M Khalat- 
nikov [ZlllEj] and V A Belinsky, E M Lifshitz and I M Khalatnikov [TO] ELD [□] ■ 
These papers revealed the enigmatic phenomenon of an oscillatory approach to 
the singularity which has become known also as the Mixmaster Universe [13] . 
The model of the closed homogeneous but anisotropic universe with three de- 
grees of freedom (Bianchi IX cosmological model) was used to demonstrate that 
the universe approaches the singularity in such a way that its contraction along 
two axes is accompanied by an expansion with respect to the third axis, and 
the axes change their roles according to a rather complicated law which reveals 
a chaotical behavior [TTI [T2l \U\ [T5 ] . 

The study of the dynamics of the universe in the vicinity of the cosmological 
singularity has become an explodingly developing field of modern theoretical 
and mathematical physics. First of all we would like to mention the generaliza- 
tion of the study of the oscillatory approach to the cosmological singularity in 
multidimensional cosmological models. It was noticed |16j that the approach to 
the cosmological singularity in the multidimensional (Kaluza-Klein type) cos- 
mological models has a chaotic character in the spacetimes whose dimensionality 
is not higher then ten, while in the spacetimes of higher dimensionalities a uni- 
verse after undergoing a finite number of oscillations enters into a monotonic 
Kasner-type contracting regime. 

The development of cosmological studies based on superstring models has 
revealed some new aspects of the dynamics in the vicinity of the singularity 
[18]. First, it was shown that in these models exist mechanisms of changing of 
Kasner epochs provoked not by the gravitational interactions but by the influ- 
ence of other fields present in these theories. Second, it was proved that the 
cosmological models based on six main superstring models plus D = 11 super- 
gravity model exhibit the chaotical oscillatory approach towards the singularity. 
Third, the connection between cosmological models, manifesting the oscillatory 
approach towards singularity and a special subclass of infinite-dimensional Lie 
algebras [H] - the so called hyperbolic Kac-Moody algebras was discovered (a 
comprehensive review of the corresponding mathematical tools with their ap- 
plication to the BKL studies was given in [ID])- The study of the algebraic 
structures, staying behind the chaotical approach to the cosmological singular- 
ity opens some new (though still very weakly elaborated) prospectives for the 
development of the consistent quantum gravity theory [3T] . 

Speaking about the new aspects of the oscillatory approach to the cosmolog- 
ical singularity in the multidimensional and superstring theories, one should not 
forget that already "classical" BKL behaviour for the 3+1 dimensional general 
relativity is not yet totally understood and worth further studying. Besides, I 
can try to attract attention to some philosophical aspects of this phenomenon 
which were underestimated until now. 
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The structure of the paper is the following: in Sec. 2 we briefly discuss the 
Landau theorem about the singularity which was not published in a separate 
paper and was reported in the book [22] and in the review [7]; in Sec. 3 we 
shall recall the main features of the oscillatory approach to the singularity in 
relativistic cosmology, including its chaoticity; Sec. 4 will be devoted to the 
modern development of the BKL ideas and methods, including the dynamics in 
the presence of a massless scalar held, multidimensional cosmology, superstring 
cosmology and the correspondence between chaotic cosmological dynamics and 
hyperbolic Kac-Moody algebras; in the conclusive Sec. 5 we shall express some 
thoughts about the past and the future of the universe in light of the BKL 
phenomenon. 



2 Landau theorem about the singularity 

Let us consider the synchronous reference frame with the metric 

ds 2 = dt 2 -~/ a pdx a dx f> , (1) 

where 7 Q ^ is the spatial metric. L D Landau pointed out that the determinant g 
of the metric tensor in a synchronous reference system must go to zero at some 
finite time provided some simple conditions on the equation of state are satisfied. 
To prove this statement it is convenient to write down the — component of 
the Ricci tensor as 

l dK« _ I 
2 dt 4' 

where K a p is the extrinsic curvature tensor defined as 



< = - 7*5*?. ( 2 ) 



K aP = — (3) 

and the spatial indices are raised and lowered by the spatial metric j a /3- The 
Einstein equation for is 

= - \t, (4) 
where the energy-momentum tensor is 

T- = (p+ y)UiV? - <r?p, (5) 

where p,p and itj are the energy density, the pressure and the four- velocity 
respectively. The quantity Tq — in the right-hand side of equation ^ is 

T°-±T=±(p + 3p) + (p + p)u a u a , (6) 
and is positive provided 

p + 3p > 0. (7) 

Thus, from Eq. (g]) follows that 

1 f)K a 1 

2^f + 4^-°- (8) 
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Because of the algebraic inequality 



K£K? > \{K2) 2 (9) 



we have 

dKZ 1 



-AKZY < o (io) 



dt 6 
or 

5 11 

> -. (11) 

dt K% ~ 6 y } 

If K" > at some moment of time, then if t decreases, the quantity decreases 
to zero within a finite time. Hence K" goes to +oo and that means because of 
the identity 

*s-^-s h « (12) 

that the determinant g goes to zero (no faster than t 6 according to the inequality 
(fTTjO . If at the initial moment K™ <, then the same result will be obtained for 
the increasing time. The similar result was obtained by Raychaudhuri [23j for 
the case of dust- like matter and by Komar [24] . 

This result does not prove that it is inevitable that there exists a true physical 
singularity, which is one belonging to the spacetime itself and is not connected 
with the character of the chosen reference system. However, it played an im- 
portant role stimulating the discussion about the existence and generality of 
the singularities in cosmology. Notice that the condition of the energodomi- 
nance ((7]) used for the proof of the Landau theorem appears also in the proof 
of the Penrose and Hawking singularity theorem [6] . Moreover, the breakdown 
of this condition is necessary for an explanation of the phenomenon of cosmic 
acceleration. 

The Landau theorem is deeply connected with the appearance of caustics 
studied by E M Lifshitz, I M Khalatnikov and V V Sudakov [25] and discussed 
between them and L D Landau in 1961. Trying to construct geometrically 
the synchronous reference frame, one start from the three-dimensional Cauchy 
surface and design the family of geodesies orthogonal to this this surface. The 
length along these geodesies serves as the time measure. It is known that these 
geodesies intersect on some two-dimensional caustic surface. This geometry 
constructed for the empty space is valid also in the presence of dust-like matter 
(p = 0). Such matter, moving along the geodesies concentrates on caustics, 
but the growth of density cannot be unbounded because the arising pressure 
destroys the caustics 0- This question was studied by L P Grishchuk [2B]. Later 
V I Arnold, S F Shandarin and Ya B Zeldovich j27J have used the caustics for the 
explanation of the initial clustering of the dust, which though does not create 
physical singularities, is nevertheless responsible for the creation of the so called 
pancakes. These pancakes represent the initial stage of the development of the 
large scale structure of the universe. 



an empty space the caustic is a mathematical, but not a physical singularity. This 
follows simply from the fact that we can always shift its location by changing the initial 
Cauchy surface. 
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3 Oscillatory approach to the singularity in rel- 
ativist ic cosmology 

One of the first exact solutions found in the framework of general relativity 
was the Kasner solution [T7] for the Bianchi-I cosmological model representing 
gravitational field in an empty space with Euclidean metric depending on time 
according to the formula 

ds 2 = dt 2 - t 2pi dx 2 - t 2p2 dy 2 - t 2 ^dz 2 , (13) 

where the exponents p\ , P2 and p^ satisfy the relations 

Pi+P2+P3=pI+pI+P 2 3 = 1- (14) 

Remarkably, this solution was the first non-stationary cosmological solution, 
found before isotropic Friedmann solutions. Perhaps, because of its "exoticity" 
it was for many years ignored by the working cosmologists and has become 
appreciated only in the fifties. 

Choosing the ordering of exponents as 

Pi <P2 < P3 (15) 

one can parameterize them as [7] 

—u l + u u(l + u) 

1 + u + u A 1 + u + u z 1 + u + u A 

As the parameter u varies in the range u > 1, pi,P2 and p% assume all their 
permissible values 

~g <Pl 0< P 2 < ^ | <P3 < I- (17) 

The values u < 1 lead to the same range of values of p\, p2, P3 since 

Pl (u) =Pl<yU ^' P2 (^) =P3 ( U )' P3 f^J = Pz( u )- ( 18 ) 

The parameter u introduced in early sixties has appeared very useful and its 
properties attract attention of researchers in different contests. For example, 
in the recent paper [55] a connection was established between the Lifshitz- 
Khalatnikov parameter u and the invariants, arising in the context of the Petrov's 
classification of the Einstein spaces [33] . 

In the case of Bianchi-VIII or Bianchi-IX cosmological models the Kasner 
regime ([IB" ]) . (TH)) ceased to be an exact solution of Einstein equations, however 
one can design the generalized Kasner solutions [5J-[I1]- It is possible to con- 
struct some kind of perturbation theory where exact Kasner solution ([IB"]) . (TH)) 
plays role of zero-order approximation while role of perturbations play those 
terms in Einstein equations which depend on spatial curvature tensors (appar- 
ently, such terms are absent in Bianchi-I cosmology). This theory of perturba- 
tions is effective in the vicinity of singularity or, in other terms, at t — > 0. The 
remarkable feature of these perturbations consists in the fact that they imply 
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the transition from the Kasner regime with one set of parameters to the Kasner 
regime with another one. 

The metric of the generalized Kasner solution in a synchronous reference 
system can be written in the form 

ds 2 = dt 2 — {a 2 l a lfj + b 2 m a mp + c 2 n a np)dx a dx l3 , (19) 

where 

a = t Pl , b = t Pm , c = t Pn . (20) 

The three-dimensional vectors I, m, n define the directions along which the spa- 
tial distances vary with time according to the power laws (1201) . Let pi = pi, p m = 
P2,Pn = P3 so that 

a~t Pl , b~t p \ c~t P3 , (21) 

i.e. the Universe is contracting in directions given by vectors m and n and is 
expanding along I. It was shown [11] that the perturbations caused by spatial 
curvature terms make the variables a, b and c to undergo transition to another 
Kasner regime characterized by the following formulae: 

a~ t p >, b~t p '\ c~t p '*, (22) 

where 

. _ bil w _ 2|pi|-p 2 , _ p 3 -2|pi| , . 

Vl ~ l-2\ Vx \> Vm - l-2bx|' P "- 1-2| P1 |- (26) 

Thus, the effect of the perturbation is to replace one "Kasner epoch" by 
another so that the negative power of t is transformed from the I to the m 
direction. During the transition the function a(t) reaches a maximum and b(t) 
a minimum. Hence, the previously decreasing quantity b now increases, the 
quantity a decreases and c(t) remains a decreasing function. The previously 
increasing perturbation caused the transition from regime ([21 p to that (|2"2"|) is 
damped and eventually vanishes. Then other perturbation begins grow which 
leads to a new replacement of one Kasner epoch by another, etc. 

We would like to emphasize that namely the fact that perturbation implies 
such a change of dynamics which extinguishes it gives us an opportunity to use 
perturbation theory so successfully. Let us note that the effect of changing of the 
Kasner regime exists already in the cosmological models more simple than those 
of Bianchi IX and Bianchi VIII. As a matter of fact in the Bianchi II universe 
there exists only one type of perturbations, connected with the spatial curvature 
and this perturbation makes one change of Kasner regime (one bounce). This 
fact was known to E M Lifshitz and I M Khalatnikov at the beginning of sixties 
and they have discussed this topic with L D Landau (just before the tragic 
accident) who has appreciated it highly. The results describing the dynamics 
of the Bianchi IX model were reported by I M Khalatnikov in his talk given 
in January 1968 in Henri Poincare Seminar, in Paris. John A Wheeler who 
was present there pointed out that the dynamics of the Bianchi IX universe 
represents a non-trivial example of the chaotic dynamical system. Later Kip 
Thorn has distributed the preprint with the text of this talk. 

Coming back to the rules governing the bouncing of the negative power of 
time from one direction to another it is necessary to emphasize that the very 
complicated system of non-linear differential equations in partial derivatives 
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is reduced in the vicinity of singularity to rather a simple system of ordinary 
differential equations. To extract the information about the rules (|23p it was 
enough to make a qualitative analysis of the latter. This analysis could be 
compared with the description of the motion of the ball climbing to the top 
of a hill: after the reaching of the highest possible point, it stops and begins 
rolling down. At the foot of the hill its velocity is equal to its initial velocity, 
but with the opposite sign. Besides, some kind of the law of the conservation of 
the sum of velocities corresponding to the expansion (contraction) of different 
space directions was used [TTJ [25] 

On the other hand, it was shown that the rules of bouncing (123)) could be 
conveniently expressed by means of the parametrization (1161) : 



and then 



pi =pi(u), p m =p 2 (u), p n =P3(u) (24) 



Pi =P2(u - 1), p' m = Pi(u - 1), p' n =p 3 (u - 1). (25) 

The greater of the two positive powers remains positive. 

The successive changes (1251) . accompanied by a bouncing of the negative 
power between the directions I and to, continue as long as the integral part 
of u is not exhausted, i.e. until u becomes less that one. Then, according to 
Eq. (1181) the value u < 1 transforms into u > 1, at this moment either the 
exponent pi or p m is negative and p n becomes smaller one of the two positive 
numbers (p n = P2)- The next sequence of changes will bounce the negative 
power between the directions ft and I or n and m. Let us emphasize that 
the usefulness of the Lifshitz-Khalatnikov parameter u is connected with the 
circumstance that it allows to encode rather complicated laws of transitions 
between different Kasner regimes ((23)) in such simple rules as u — > u — 1 and 

u 

Consequently, the evolution of our model towards a singular point consists of 
successive periods (called eras) in which distances along two axes oscillate and 
along the third axis decrease monotonically, the volume decreases according to 
a law which is near to ~ t. In the transition from one era to another, the 
axes along which the distances decrease monotonically are interchanged. The 
order in which the pairs of axes are interchanged and the order in which eras of 
different lengths follow each other acquire a stochastic character. 

To every (sth) era corresponds a decreasing sequence of values of the param- 
eter u. This sequence has the form Umax, Umax - 1, ... , u m \ n , where u^ in < 1. 
Let us introduce the following notation: 

<41=z (s) , u[± x = k^+x^ (26) 



i.e. fcW = [«mL] (the square brackets denote the greatest integer < 
The number k^ defines the era length. For the next era we obtain 

(s+l) = _L t,(s+l) = 
max ( c\ i n 



■ma:r } 



U 

X 



max 



5 W 



(27) 



The ordering with respect to the length of of the successive eras (mea- 
sured by the number of Kasner epochs contained in them) acquires asymptoti- 
cally a stochastic character . The random nature of this process arises because 
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of the rules (|2l))) -(f2~T f which define the transitions from one era to another in 
the infinite sequence of values of u. If all this infinite sequence begins since 
some initial value Umax = k^ + x^°\ then the lengths of series k^°\ k^ l \ . . . are 
numbers included into an expansion of a continuous fraction: 

fc(°) + x «» = fc(°) + - 1 , (28) 

k {1 > + jmr^ 

We can describe statistically this sequence of eras if we consider instead of a 
given initial value Umax = k^ +x^°' > a distribution of x ^ over the interval (0, 1) 
governed by some probability law |14j . Then we also obtain some distributions of 
the values of x^ which terminate every sth series of numbers. It can be shown 
that with increasing s, these distributions tend to a stationary (independent 
of s) probability distribution w(x) in which the initial value x^ s ' is completely 
"forgotten" : 

W W = noAi o - (29) 
(1 + x) m2 

It follows from Eq. (|29l) that the probability distribution of the lengths of series 
k is given by 

w <*> (30) 

The source of stochasticity arising at the oscillatory approach to the cos- 
mological singularity can be described in such terms: the transition from one 
Kasner era to another is described by the transformation of the interval [0, 1] 
into itself by the formula 

Tx=(-\, i.e., Xs+1 = \ — }. (31) 



x J v x , 

This transformation is expanding and possesses the property of exponential 
instability. It is not one to one transformation. Its inverse is not unique. In 
other words, fixing the value of the parameter u we can predict the evolution 
towards singularity, but we cannot describe the past. 
One can pass over from a one-sided infinite sequence 

(x ,x 1 ,x 2 ,...), (32) 

to a doubly infinite sequence [15 

X = (.. . ,X-2,X-i,XQ,X\,X2, ■ ■ ■)■ (33) 
The sequence X is equivalent to a sequence of integers 

K = (...,fc_ 2 ,&-i,ft D ,fci,Aa,...) J (34) 



such that 



Inversely, 



1 



X s -1 



1 



ks+i + 



(35) 
(36) 
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x+ = [k.,k a+1 ,...] (37) 

and 

x~ = [fc s _i, fc s _ 2 , ■ • ■] (38) 
The shift of entire sequence X to the right means a joint transformation 



+1 ~~ \xf J ' ~~ 

This is a one-to-one mapping in the unit square, which permits to calculate 
in an exact manner probability distributions for other parameters describing 
successive eras such as parameter 6 giving relation between the amplitudes of 
logarithms of functions a,b,c and logarithmic time |15) . 

Thus, we have seen from the results of statistical analysis of evolution in the 
neighbourhood of singularity [14] that the stochasticity and probability distri- 
butions of parameters arise already in classical general relativity. 

At the end of this section a historical remark is in order. The continuous 
fraction ((28)) was shown in 1968 to I M Lifshitz (L D Landau has already left) 
and he immediately noticed that one can derive the formula for a stationary 
distribution of the value of x (f2"9"|) . Later it becomes known that this formula 
was derived in nineteenth century by Karl F Gauss, who had not published it 
but had described it in a letter to one of colleagues. 



1 


+ x s ) 







(39) 



4 Oscillatory approach to the singularity: mod- 
ern development 

The oscillatory approach to the cosmological singularity described in the preced- 
ing section was developed for the case of an empty spacetime. It is not difficult 
to understand that if one considers the universe filled with a perfect fluid with 
the equation of state p — wp, where p is the pressure, p is the energy density 
and w < 1, then the presence of this matter cannot change the dynamics in the 
vicinity of the singularity. Indeed, using the energy conservation equation one 
can show that 

= Po = Po (40) 
9 ~ (abc)^ 1 ~ t w + x ' 1 ' 

where po is a positive constant. Thus, the term representing the matter in 
the Einstein equations behaves like ~ \/t 1+w and at t — > is weaker than 
the terms of geometrical origin coming from the time derivatives of the metric, 
which behaves like l/t 2 , let alone the perturbations due to the presence of 
spatial curvature, responsible for changes of a Kasner regime,which behave like 
l/£ 2+4 l pi L However, the situation changes drastically, if the parameter w is 
equal to one, i.e. the pressure is equal to the energy density. Such kind of 
matter is called "stiff matter" and can be represented by a massless scalar field. 
In this case p ~ l/t 2 and the contribution of matter is of the same order as 
main terms of geometrical origin. Hence, it is necessary to find a Kasner type 
solution, taking into account the presence of terms, connected with the presence 
of the stiff matter (a massless scalar field) . Such study was carried out in paper 
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[50] . It was shown that again the scale factors a, b and c can be represented as 
t 2pi , t 2p2 and t 2p3 respectively, where the Kasner indices satisfy the relations: 



Pi +P2 +P3 = 1, Pi 



p\ 



2 -, 2 

Pi = 1 - 9 ) 



(41) 



where the number q 2 reflects the presence of the stiff matter and is bounded by 



Q 2 < 



(42) 



One can see that if q 2 > 0, then exist combinations of the positive Kasner 
indices, satisfying the relations (|4"Tj) . Moreover, if q 2 > | only sets of three 
positive Kasner indices can satisfy the relations (|4"Tj) . If a universe finds itself 
in a Kasner regime with three positive indices, the perturbative terms, existing 
due to the spatial curvatures are too week to change this Kasner regime, and 
thus, it becomes stable. That means, that in the presence of the stiff matter, 
the universe after a finite number of changes of Kasner regimes finds itself in 
a stable regime and oscillations stop. Thus, the massless scalar field plays 
"anti-chaotizing" role in the process of the cosmological evolution [3D] . One can 
use the Lifshitz-Khalatnikov parameter also in this case. The Kasner indices 
satisfying the relations (|4l"j) are conveniently represented as [30] 



Pi = 

P2 = 

P3 = 
P 2 = 



1 + u + u 2 ' 

1 + u 
1 + u + u 2 

1 + U 14 

1 + u + u 2 

2(1 + u + u 2 ) 2 



u-l 
2 

u- 1 



(i- ( i-/?Y /2 ) 

(1-(1-/3 2 ) 1/2 ) 



(u 2 



1) 



(43) 



The range of the parameter u now is — 1 < u < 1, while the admissible values 
of the parameter q at some given u are 



q 2 < 



(u 2 - l) 2 



2(l + u + u 2 ) 2 ' 



(44) 



One can easily show that after one bounce the value of the parameter q 2 changes 
according to the rule 



9/2 2 

q -> q = q X 



1 



(45) 



Thus, the value of the parameter q 2 grows and, hence, the probability to find all 
the three Kasner indices to be positive increases. It confirms again the statement 
that after a finite number of bounces the universe in the presence of the massless 
scalar field finds itself in the Kasner regime with three positive indices and the 
oscillations stop. 

In the second half of eighties a series of papers was published [16], where 
were studied the solutions of Einstein equations in the vicinity of singularity 
for d + 1-dimensional spacetimes. The multidimensional analog of a Bianchi-I 
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universe was considered, where the metric is a generalized Kasner metric: 



ds 2 =dt 2 -J2t 2pi dx i2 , (46) 
i=l 

where the Kasner indices pi satisfy the conditions 

J> = X>? = 1. (47) 

i=l »=1 

In the presence of spatial curvature terms the transition from one Kasner epoch 
to another occurs and this transition is described by the following rule: the new 
Kasner exponents are equal to 



with 



Pi, p'n • • • i Pd = ordering of (q 1 ,q 2 , . . . ,q d ), (48) 



—Pi - P P2 



l + 2 Pl +P ' " l + 2 Pl + P' ' 

Pd-2 2p 1 + P+p d ^ 1 

c ld-2 — : ; — =J, — 



l + 2 Pl +P' l + 2 Pl +P 

31 + P + Pd 

1 + 2pi + P ' 



2 Pl +P + Pd 

Id = , I „ , 5 ' ( 49 ) 



where 



d-2 

P = ^ft. (50) 

i=2 

However, such a transition from one Kasner epoch to another occurs if at least 
one of the numbers a^fc is negative. These numbers are defined as 

a ijk = 2p t + ^2 PU (i^jJ^ k,j ^ k). (51) 

For the spacetimes with d < 10 one of the factors a is always negative and, 
hence one change of Kasner regime is followed by another one, implying in 
such a way the oscillatory behaviour of the universe in the neighbourhood of 
the cosmological singularity. However, for the spacetimes with d > 10 exist 
such combinations of Kasner indices, satisfying Eq. (|47[) and for which all the 
numbers aijk are positive. If a universe enters into the Kasner regime with such 
indices, ( so called "Kasner stability region") its chaotical behaviour disappears 
and this Kasner regime conserves itself. Thus, the hypothesis was forward that 
in the spacetimes with d > 10 after a finite number of oscillations the universe 
under consideration finds itself in the Kasner stability region and the oscillating 
regime is replaced by the monotonic Kasner behaviour. 

The discovery of the fact that the chaotic character of the approach to the 
cosmological singularity disappears in the spacetimes with d > 10 was unex- 
pected and looked as an accidental result of a game between real numbers satis- 
fying the generalized Kasner relations (|49[) . Later it becomes clear that behind 
this fact there is a deep mathematical structure, namely, the hyperbolic Kac- 
Moody algebras. Indeed, in the series of works by Damour, Henneaux, Nicolai 
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and some other authors (see e.g. Refs. [TB]) on the cosmological dynamics in 
the models based on superstring theories, living in 10-dimensional spacetime 
and on the d + 1 = 11 supergravity model, it was shown that in the vicinity of 
the singularity these models reveal oscillating behaviour of the BKL type. The 
important new feature of the dynamics in these models is the role played by 
non-gravitational bosonic fields (p-forms) which are also responsible for transi- 
tions from one Kasner regime to another. For description of these transitions 
the Hamiltonian formalism |13) becomes very convenient. In the framework of 
such formalism the configuration space of the Kasner parameters describing the 
dynamics of the universe could be treated as a billiard while the curvature terms 
in Einstein theory and also p-form's potentials in superstring theories play role 
of walls of these billiards. The transition from one Kasner epoch to another is 
the reflection from one of the walls. Thus, there is a correspondence between 
rather complicated dynamics of a universe in the vicinity of the cosmological 
singularity and the motion of an imaginary ball on the billiard table. 

However, there exists a more striking and unexpected correspondence be- 
tween the chaotical behaviour of the universe in the vicinity of the singularity 
and such an abstract mathematical object as the hyperbolic Kac-Moody alge- 
bras [T5] . Let us explain briefly what does it mean. Every Lie algebra is defined 
by its generators fti, ej, fi,i = 1, . . . ,r, where r is the rank of the Lie algebra, i.e. 
the maximal number of its generators hi which commutes each other (these gen- 
erators constitute the Cartan subalgebra) . The commutation relations between 
generators are 

[e,, fj] = Sijhi, 
[hi,ej] = Aij = AijBj, 

[hit fj] — ~Aijfj, 

[h i ,h j ]=0. (52) 

The coefficients Aij constitute the generalized Cartan r x r matrix such that 
An — 2, its off-diagonal elements are non-positive integers and Aij = for 
i =/= j implies Aji = 0. One can say that the ej are rising operators, similar to 
well-known operator L + = L x + iL y in the theory of angular momentum, while 
fi are lowering operators like L_ = L x — iL y . The generators hi of the Cartan 
subalgebra could be compared with the operator L z . The generators should 
also obey the Serre's relations 

(ad eiY'^e-j = 0, 

(adf^-^f, = 0, (53) 

where (adA)B = [A, B]. 

The Lie algebras Q{A) build on a symmetrizable Cartan matrix A have been 
classified according to properties of their eigenvalues: 
if A is positive definite, G(A) is a finite-dimensional Lie algebra; 
if A admits one null eigenvalue and the others are all strictly positive Q(A) is 
an Affine Kac-Moody algebra; 

if A admits one negative eigenvalue and all the others are strictly positive, Q{A) 
is a Lorentz KM algebra. 

There exists a correspondence between the structure of a Lie algebra and a 
certain system of vectors in the r-dimensional Euclidean space, which simplify 



12 



essentially the task of classification of the Lie algebras. These vectors called 
roots represent the rising and lowering operators of the Lie algebra. The vectors 
corresponding to the generators e\ and fi are called simple roots. The system 
of simple positive roots (i.e. the roots, corresponding to the rising generators 
ej) can be represented by nodes of their Dynkin diagrams, while the edges 
connecting (or non connecting) the nodes give an information about the angles 
between simple positive root vectors. 

An important subclass of Lorentz KM algebras can be defined as follows: A 
KM algebra such that the deletion of one node from its Dynkin diagram gives 
a sum of finite or affine algebras is called an hyperbolic KM algebra. These 
algebras are all known. In particular, there exists no hyperbolic algebra with 
rank higher than 10. 

Let us recall some more definitions from the theory of Lie algebras. The 
reflections with respect to hyperplanes orthogonal to simple roots leave the 
systems of roots invariant. The corresponding finite-dimensional group is called 
Weyl group. Finally, the hyperplanes mentioned above divide the r-dimensional 
Euclidean space into regions called Weyl chambers. The Weyl group transform 
one Weyl chamber into another. 

Now, we can briefly formulate the results of approach [T8] following the pa- 
per [in] : the links between the billiards describing the evolution of the universe 
in the neighbourhood of singularity and its corresponding Kac-Moody algebra 
can be described as follows: 

the Kasner indices describing the "free" motion of the universe between the 
reflections from the wall correspond to the elements of the Cartan subalgebra 
of the KM algebra; 

the dominant walls, i.e. the terms in the equations of motion responsible for the 
transition from one Kasner epoch to another, correspond to the simple roots of 
the KM algebra; 

the group of reflections in the cosmological billiard is the Weyl group of the KM 
algebra; 

the billiard table can be identified with the Weyl chamber of the KM algebra. 

One can imagine two types of billiard tables: infinite such where the linear 
motion without collisions with walls is possible (non-chaotic regime) and those 
where reflections from walls are inevitable and the regime can be only chaotic. 
Remarkably, the Weyl chambers of the hyperbolic KM algebras are designed in 
such a way that infinite repeating collisions with walls occur. It was shown that 
all the theories with the oscillating approach to the singularity such as Einstein 
theory in dimensions d < 10 and superstring cosmological models correspond to 
hyperbolic KM algebras. 

The existence of links between the BKL approach to the singularities and the 
structure of some infinite-dimensional Lie algebras has inspired some authors 
to declare a new program of development of quantum gravity and cosmology 
pH] , They propose "to take seriously the idea that near the singularity (i.e. 
when the curvature gets larger than the Planck scale) the description of a spa- 
tial continuum and space-time based (quantum) field theory breaks down, and 
should be replaced by a much more abstract Lie algebraic description. Thereby 
the information previously encoded in the spatial variation of the geometry and 
of the matter fields gets transferred to an infinite tower of Lie- algebraic vari- 
ables depending only on "time" . in other words we are led to the conclusion 
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that space-and thus, upon quantization also space-time - actually disappears 
(or "de-emerges") as the singularity is approached." 

5 Conclusion: some thoughts about the past 
and the future of the universe 

In the preceding section we have told briefly about the newest development of 
the theory of the BKL approach to the cosmological singularity connected with 
the superstring-inspired cosmological models and the infinite-dimensional Lie 
algebras. However, already in the "standard" 3+1 - dimensional general rela- 
tivity the effect of the oscillatory approach to the singularity and the chaoticity 
implied by it is of great interest. Indeed, the discovery of the non-static time- 
dependent cosmological solutions in the general relativity, first of all the Fried- 
mann solutions, has given birth to the animate discussions on such questions as: 

Did the universe have the beginning ? 
Will the universe have the end ? 

Can the Universe exist during a finite interval of time ? 

What was before the beginning and what will be after the end ? 

These questions look quite reasonable, because we know that in all three Fried- 
mann models -flat, open and closed - the universe has the beginning and this be- 
ginning is nothing but the Big Bang singularity. In the closed Friedmann model 
the universe has also the end - the Big Crunch singularity and exists during a 
finite period of time. Moreover, according to the so called Standard Cosmologi- 
cal Model, based on a rather large set of observational data, something like Big 
Bang has taken place approximately 13.7 billions years ago (measured in terms 
of the cosmic, i.e. synchronous time). The more or less accepted existence of the 
beginning of the evolution of the universe and the possible existence of the end 
of the universe can be source of the joy for those who believe into the creation 
of the universe and for whom its possible end can also be confirmation of their 
philosophical or theological beliefs. It is curious that the Pontifical Academy of 
Sciences has organized a special conference in the October-November of 2008 in 
Vatican with the title " Scientific insights into the evolution of the universe and 
of life" . On the other hand, the possibility of the finite-time existence of the 
universe can provoke some kind of psychological discomfort in persons whom 
the finite duration of the existence of the universe seems to be senseless. For 
somebody, the fact that its own existence takes place in the universe which lives 
only during a finite period of time, can look depressive. 

Analyzing this aspects of the problem of the evolution of the universe, one 
should ask himself /herself, which time parametrization should we use, speaking 
about the time of the existence of the Universe ? As we know since the epoch 
of the creation of the special relativity, the time is relative. In the framework of 
the general relativity the time becomes even more relativity and can run with 
different rates at different spatial points. Making conformal transformations 
(for example, then one constructs the Penrose conformal diagrams [B]) one can 
turn an infinite time interval into the finite time interval. Thus, why should 
we use the cosmic time ? The answer on this question is simple: the cosmic 
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time for a particle staying in rest in a Friedmann homogeneous and isotropic 
universe coincides with the proper time, introduced in the special relativity. 
Hence, when we are considering the present day universe it is quite reasonable 
to discuss it in terms of the cosmic time and to say that the universe was created 
13.7 billions years ago. However, when we consider the vicinity of the Big Bang 
cosmological singularity in the past, or when we accept the possibility of exis- 
tence of the future Big Crunch singularity in the distant future, the situation 
changes drastically. The universe is extremally anisotropic in the neighbour- 
hood of such singularities and is described by the chaotic succession of Kasner 
epochs and eras, as was discussed above. (One can make a precision here: while 
choosing very special isotropic initial conditions one can avoid the chaoticity in 
the neighbourhood of the Big Bang singularity, which can , in principle have 
the Friedmannian form, it is impossible not to have the chaotic regime in the 
vicinity of the Big Crunch singularity, because, the inhomogeneities, developed 
during the evolution of the universe, make its contracting stage highly anistropic 

BZD- 

Thus, while the evolution between an arbitrary moment of cosmic time and 
the moment corresponding to the initial Big Bang or final Big Crunch singularity 
occupies a finite interval of cosmic time, an infinite number of events occurs 
during this finite period. The infinite chaotic succession of Kasner epochs and 
eras, makes senseless the cosmic time as a measure of a cosmological evolution. 
Indeed, we have an infinite history that separates us from the birth of the 
Universe at the Big Bang. If the contraction of the universe culminating in the 
encounter with the Big Crunch singularity waits for us in the future [53], we 
still have an infinite number of events in front of us. Thus, the BKL oscillatory 
regime of approaching to the cosmological singularity screens us from the Big 
Bang and Big Crunch. 

From the mathematical point of view it means, that in the vicinity of sin- 
gularity the natural time parameter is not the cosmic time, but the logarithmic 
time. Then, while the cosmic time runs from the moment zero, corresponding 
to the singularity to some finite moment t\, the logarithmic time runs from — oo 
to lnti, spanning an infinite interval of time. 

Remarkably, the comment concerning the importance of the logarithmic time 
can be found already in the last before last paragraph of the Landau and Lif- 
shitz monograph [35]: "The successive series of oscillations crowd together as 
we approach the singularity. An infinite number of oscillations are contained 
between any finite world time t and the moment t = 0. The natural variable 
for describing the time behaviour of this evolution is not the time t, but its 
logarithm, In t, in terms of which the whole process of approach to the singular 
point is spread out to -co". 

The similar idea is expressed also in the cited earlier paper by Damour and 
Nicolai [2T]: "There is no "quantum bounce" bridging the gap between an in- 
coming collapsing and an outgoing expanding quasi-classical universe. Instead 
"life continues" at the singularity for an infinite affine time,but with the under- 
standing that (i) dynamics no longer "takes place" in space, and (ii) the infinite 
affine time [measured, say, by the Zeno-like time coordinate t] corresponds to a 
sub-Planckian interval < T < Tpi anc k of geometrical proper time." Curiously, 
the analog of the object that the authors of [2T] call Zeno-like time is the so 
called spatial tortoise coordinate in the Schwarzschild geometry [53] . Both these 
names have their origin in the Zeno's paradox about Achilles and the tortoise, 
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which is, perhaps, the first example of the transforming of the finite time in- 
terval into an infinite one (see, e.g. the section I of the third part of the third 
volume of the "War and Peace" by Leo Tolstoj [55]). 

Concluding, I would like to say semiseriously that the discovery of the os- 
cillatory approach to the cosmological singularity has a practical meaning: it 
liberates us from the fear of the End of the World. 

This work was partially supported by the RFBR grant 08-02-00923 and by 
the grant LSS-4899.2008.2. 
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